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Abstract 

1^ ' We formulate general conditions on 3 x 3 neutrino mass matri- 

ces under which a degenerate pair of neutrinos at a high scale would 
^ , split at low scale by radiative corrections involving only the standard 

\ model fields. This generalizes the original observations of Wolfenstein 

on pseudo Dirac neutrinos to three generations. A specific model in- 
volving partially broken discrete symmetry and solving the solar and 
atmospheric anomalies is proposed. The symmetry pattern of the 
model naturally generates two large angles one of which can account 
for the large angle MSW solution to the solar neutrino problem. 

1 Introduction 

Dirac neutrinos are associated with an unbroken U(l) symmetry acting on 
leptons. A small breaking of this symmetry splits a Dirac neutrino into a 
pair of majorana neutrinos with (mass)^ difference much smaller than the 
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square of the original mass. Such a pair can simultaneously describe small 
splitting and large mass and is of phenomenological importance in solving 
solar and/or atmospheric neutrino anomalies 

Zeroth order approximation to a pseudo Dirac neutrino is provided by 
the following texture in case of two generations say, e and fi 



This displays an unbroken Le—L^ symmetry and its breaking ( by introducing 
small non-zero diagonal elements) leads to a pseudo Dirac neutrino. This 
breaking can be explicitly introduced by allowing for additional fields like 
Higgs triplet or right handed neutrinos or it can be introduced radiatively 
by breaking the symmetry Lg — in the charged lepton sector Q . 

A slightly non-trivial example of the pseudo Dirac neutrino is provided 
by a mass matrix discussed by Wolfenstein H 





(2) 



Both the textures in eqs.(l) and (2) lead to a pair of neutrinos with equal 
and opposite eigenvalues. However they differ conceptually and phenomeno- 
logically from each other Q Unlike in eq.(l), the matrix in eq.(2) cannot 
be invariant under a U{1) symmetry corresponding to any combination of 
lepton numbers. This has the consequence that the charged current interac- 
tions defined in the mass basis of the degenerate pairs violate lepton number 
[|] and the pseudo Dirac pair gets split automatically by radiative correc- 
tions 0. Thus the theory described by eq. (2) intrinsically defines a pseudo 
Dirac neutrino while one needs to invoke additional fields in order to break 
the Lg — symmetry in case of eq.(l). At the phenomenological level, the 
mixing implied by (1) is maximal while it is arbitrary in case (2). The most 
phenomenological discussions related to pseudo Dirac neutrinos in the liter- 
ature are in the context of texture in (1). We wish to discuss here instead 
several interesting aspects related to the Wolfenstein texture, eq.(2) and its 
generalization to three families. 

^ We are assuming here that these textures are defined in a basis with a diagonal 
charged lepton mass matrix. 
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The radiative splitting of neutrinos in case (2) is most simply demon- 
strated using the relevant renormalization group (RG) equations. Assume 
that the neutrino mass matrix in is specified at some high scale Mx and 
the effective theory below this scale is the standard model (SM) or the mini- 
mal supersymmetric standard model (MSSM). The neutrino mass matrix at 
a low scale /i ~ Mz is then given by [0 

MM ^ PM,{Mx)P , (3) 

where in the three-generation case P = Diag.{l + 6e,l + 5^, 1 + Sr) and 

S = In^ (4) 

in case of the minimal supersymmetric standard model. ma,v ~ 174 GeV 
here refer to the charged lepton masses and the weak scale respectively. As 
before, M,^ is specified in the physical basis of the charged leptons. The 
texture (2) at high scale gets transformed to the texture 



a(l + 25,) 6(1 + 4 + 5^) 
6(l+5e + M -a(l + 25^) 



+ 0(6') (5) 



at the low scale. This describes a split pair of neutrinos. In contrast, eq.(l) 
leads to a degenerate pair even after RG evolution as in eq.@ is taken into 
account. 

The modified mixing angle and the mass splitting implied by eq.(|^) are 
given by 

tan2^(;u) ^ tan2e{Mx){l + 0{6^)) , 

A ^ 4m^cos2%. (6) 



with mo = and tan26' = ^. As follows from eq.(§), typical strength 

of radiative corrections is 5^ ~ 10^^. From phenomenological point of view, 
the required value of A can be as small as IQ-^^eVl If one starts with a 
mass matrix having such splitting at Mx then there is a possibility that the 
radiative corrections may lead to much larger splitting than this. Likewise, 
maximal mixing angle at Mx can also get destabilized 0. In the present 
case, splitting is zero at Mx and it is only induced by radiative corrections. 
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As follows from eq.(^, the magnitude of this splitting can be in the range 
10~^° — 10~^^ eV^ for mo near the atmospheric neutrino scale. Moreover, the 
mixing angle at Mx is arbitrary and receives corrections only at 0(5^) in this 
case. Thus, the texture in (2) is stable against radiative corrections unlike 
some of the textures discussed in . These properties make the Wolfenstein 
texture in (Q) also phenomenologically realistic. Theoretically, this texture 
is not a fine tuned possibility but can arise from imposition of the following 
discrete symmetry on neutrino mass matrix: 

Le iL^ ; — >■ —iLe (7) 

This symmetry is broken by hierarchical charged lepton masses which lead 
to radiative splitting between the degenerate pair. 

The 2x2 texture of eq.(^ is successful but not complete from the point of 
view of simultaneous solution to the solar and atmospheric neutrino anoma- 
lies. This would require going beyond two generations. The purpose of this 
note is to generalize the above considerations to the realistic case of three 
generations and identify phenomenologically viable models/textures leading 
to pseudo Dirac neutrinos. We first write down the general conditions on 
an arbitrary 3x3 neutrino mass matrix under which it leads to a pair of 
degenerate neutrinos. Such matrices can be classified in two categories, those 
in which degeneracy is preserved by radiative corrections involving standard 
model fields and those in which the theory describes a pseudo Dirac state. 
We formulate general criteria to distinguish between these two categories 
and show that they can be identified by looking at the structure of the lep- 
tonic mixing matrices implied at the tree level. Then we discuss an example 
which satisfies phenomenological requirements to obtain a solution to solar 
and atmospheric neutrino anomalies. Starting with a generalization of the 
original Wolfenstein mass matrix, splitting of neutrino states needed for solar 
neutrino anomaly arises in this example through radiative corrections. The 
model can lead to either vacuum solution with bi-maximal mixing or MSW 
solution corresponding to large mixing angle. 
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2 Pseudo Dirac neutrinos: General analysis 

Let us consider a CP conserving theory specified by a general 3x3 real 
symmetric mass matrix Mi, for the neutrinos: 

1 



-C„, = -{u'aLY{M,)a(3iy'pL + H.C. . (8) 

The Ml, would contain two equal and opposite eigenvalues if it satisfies the 
following condition: 

tr{M,)Y,\ = detM,, (9) 

i 

where Aj represents the determinant of the 2x2 block of M^, obtained 
by blocking i*^ (i=l,2,3) row and column. Only M^, satisfying condition 
@ would lead to a Dirac or pseudo Dirac neutrinos. Define a which 
diagonalizes such M^: 

W MM'"^ = Diag.{m, -m, m') . (10) 

The physical mass basis for neutrinos is defined by i>i = W^u'i. The neutrino 
masses can be written in terms of Majorana spinors 

t^l,3 = l^l,3L + {l^l,3Ly ^2 = V2L - {^2lY ■ (H) 

Explicitly, 

- Cm = - [m{uiui + Z/2Z/2) + m'z/sz/s] . (12) 
The degeneracy of two of the Majorana states allows us to define 

^ = ^K + ^2) (13) 

and rewrite eq.(ll2|) as 



- Cm = mipip + ^m'usus . (14) 

Note that the V' is a four component Dirac field since ip'^ i^- As a result, the 
system specified by eq.(^ corresponds to a Dirac and a majorana neutrino. 
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The charged current interactions can be written as follows in the leptonic 
mass basis: 



Here CaL represents the physical mass basis for the charged leptons and the 
K represents the leptonic Kobayashi Maskawa (KM) matrix. 

Eq.(|T^) is a straightforward generalization of the 2x2 case considered 
in 1^. It shows that although mass term for ip is invariant under a U(l) 
symmetry the charged current violates it and the Dirac state will split by 
radiative corrections 0. Thus any 3x3 matrix satisfying condition @ at 
tree level would generically lead to a pseudo Dirac state. 

While the lepton number violation is generically present, it is easy to 
identify all the 3x3 structures for M^^ which admit an unbroken f/(l) sym- 
metry corresponding to a truly Dirac neutrino. Necessary condition for this 
to happen is easy to write down using eg . ([I5|) : 

K^i = ri^K^2 (16) 

Here 7]^ = ±1. The above equation ensures that either only ip or ip'^ couples to 
a given charged lepton e^. As a result, phase rotation of ip can be symmetry 
of eq.(|l5D. 

The above condition is also sufficient to ensure truly Dirac state in case 
of two generations. But the couplings of the third generation to the first 
two requires additional constraint to obtain a U{1) symmetry. To see this, 
note that orthogonality of K does not allow all the rja in eq.(^) to have the 
same sign. Without loss of generality we can choose rji = —ri2 = —1]^ = 1. 
Other solutions of eq.(|T^) are obtained by interchange 1 2, 1 ^ 3 or an 
overall multiplication of rja by -1 in all these cases. It is easy to write a 
parameterization of K with this choice of 77 using orthogonality. 

/ ^ ^ \ 

K = ^ e 

\ V2 V2 ^ J 



(17) 
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Apart from interchange of rows and overall multiplications of any column 



by -1, this is the most general form for any real K satisfying eq.(16). Using 
above eq. ( [T7| ) in (|I^) we immediately see that if s = 0, the third neutrino 
does not mix with the first two and eq. (|15D is invariant under the symmetry, 

In this case, the third neutrino is decoupled and we get an unbroken — 
symmetry and a truly Dirac neutrino. 

Even when s ^ 0, W interactions can still be made invariant under the 
following symmetry 

But this symmetry would be violated by the mass of the third neutrino and 
hence we would expect splitting of the neutrino pairs radiatively in this case 
even though the charged current interactions possess an unbroken symmetry. 

It follows from the above discussion that there are two ways in which 
a degenerate pair at high scale would be split by the radiative corrections 
involving W: (i) When mixing matrix does not satisfy the necessary condition 
in ( p!6D and (ii) when mixing matrix satisfies this condition but the third 
neutrino has a mass and is not decoupled from the first two. The latter 
situation is more interesting since in this case the implied form of K, eq.(|T7|) 
corresponds to bi-maximal mixing which is argued [|, |^ to be useful in 
simultaneous solution of the solar and atmospheric neutrino anomalies. The 
mass of the third generation plays a non-trivial role in splitting the degenerate 
pair in this case. This was found to be true in specific example considered 
in |T0[. The present discussion highlights the importance of a non-zero third 
generation mass from more general considerations. 

All the above conclusions were based on the U{1) invariance of the W 
interactions displayed in eq. (plSf ). The same conclusions can be drawn from 
the study of the RG evolution of the texture as in eq.(0). The matrix K 
in eq.([T7|) implies the following texture for the neutrino mass matrix in the 
charged lepton mass basis: 



/ —mc ms \ 

—mc m's^ m'cs 
\ ms m'cs m'(? ) 



(18) 
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Assuming this texture to be true at Mx-, texture at a lower scale can be 
worked out using eq.(^ 



-mc(l + 5e + 5^) 
\ ms{l + 5e + 5r) 



—mc(l + Se + S^j 



m s 



-m'cs(l + Sr + 



ms(l + Se + Sr) \ 

-m'cs(l + 6r + 6fj,) 
m'c\l + 26r) J 



(19) 



This texture is seen to lead to a degenerate pair at low scale if s = or m' = 
in accordance with the conditions discussed above. Otherwise, it describes 
a pseudo Dirac state. In the former case, the neutrino mass matrix in the 
charged lepton mass basis respects either Lg — Lf^ (s = 0) or Le — — Lr 
[m' = 0) symmetry. In the latter case, splitting occurs but only at 0(5^) 
as can be seen using eq.([T9|). This makes the case (ii) discussed above also 
more natural from the point of view of stability as argued in [ITU . 



3 A Specific Model 

A very economical scheme for understanding the solar and atmospheric neu- 
trino anomaly is provided by the following scenario The neutrino 
spectrum at tree level consists of one massless neutrino and two degenerate 
neutrinos with mass in the atmospheric neutrino range. The radiative cor- 
rections split this pair and provides the scale needed to understand the solar 
neutrino anomaly. Different possibilities realizing this scenario have been 
proposed [0, We give here a specific and very economical example in the 
context of pseudo Dirac texture where one does not need to invoke any new 
physics and the W interactions provide a source for the solar scale. 

As we indicated in introduction, the Wolfenstein type structure as in 
eq.(0) can result from a discrete symmetry. A similar symmetry can be used 
to obtain a 3 x 3 generalization of eq. (0). Consider, 

{L'e, L'^) i{L'^, L'J ; (L^, e'^, fi'^, r'j^) -i{L'^, fi'n, e'n, t'j^) (20) 

Here, L'^ and {e'j^, /x'^, rjj) denote the leptonic doublets and singlets respec- 
tively. All other fields are assumed neutral with respect to the above sym- 
metry. The neutrino mass matrix invariant under this symmetry has the 
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following form: 

fa b\ 
M^{Mx) = -a b (21) 
V 6 b ) 

This mass matrix satisfies condition, eq.(^ and hence leads to a degenerate 
pair of neutrinos with mass m = V ct^ + 26^ at the tree level. It is not in 
the basis with diagonal charged leptons. The symmetry in eq.(pi]|) allows the 
following general form for the charged lepton masses: 



M, 



/ mi 1712 
—1712 —rrii 

\ 777,4 "^4 



7713 \ 
-7723 
7775 / 



(22) 



Here 777j are parameters arising in the standard way through Yukawa cou- 
plings of leptons to the Higgs field. The Mi can be diagonalized in the 
standard manner. Eg . (^21) leads to 



MiMj 



(A 
B 



B 
A 

-C 



C \ 

-c 

D I 



The parameters in the above matrix can be read off from eq.(|2 



Diag.{ml,ml,ml] 



With 



= i?23(0)i?i2(vr/4) 



1 

\ V2 



1 

V2 



\ 

C0 J 



(23) 

One finds, 
(24) 

(25) 



where C(f, = coscf)] = sine/) and 



tan ! 



2^/2C 



D-A + B 

While fairly large ranges can be allowed for the parameters 777j, we will make 
the following choice for illustrative purpose: 



nil ~ "^2 



7773 ~ 0(777^) ; 7775 ~ 0(777^) 



(26) 
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It can be seen that this choice reproduces the charged lepton masses correctly. 
Moreover, with this choice, 



tan 20 ~ 0{ 



m. 



(27) 



The neutrino mass matrix in eq. (|2l|) is diagonahzed by 



cos^l 



sin^ I 



sm 
1 

'V2 



2 

sin 6* 



cos 



2 e 



1 

V2 



2 

sin 6 



^sine \ 
-^sin^^ 



cos^ 



(28) 



with tan^ = 

a 

In spite of the complicated form for U^, the KM matrix K = U^U^'^ has 
a bi- maximal form given in eq. (|r^) . This form coupled with the masslessness 
of the third generation ensures that the degenerate pair does not split radia- 
tively. This is not surprising since the charged lepton mass matrix, eq.(P2|) is 
also invariant under the symmetry of eq.(pOD which was responsible for the 
degeneracy of neutrinos. We need to break this symmetry in order to obtain 
splitting of the degenerate pairs. One possible breaking is as follows: 



(29) 



Such a term can arise from Yukawa couplings with the standard Higgs in 
which case it breaks the symmetry explicitly. This type of hard breaking 
advocated in several papers 0] can make the theory non-renormalizable. 
This is avoided by Yukawa couplings of an additional Higgs field which is 
odd under the discrete symmetry. The discreet symmetry would then be 
spontaneously broken. Alternatively, we can work with the hard breaking in 
( |29D but assume the model to be embedded in the supersymmetric theory 
which would not make the model non-renormalizable. We shall choose the 
latter alternative in what follows. 

The inclusion of the term in eq. (|29|) , would be expected to split the degen- 
eracy between the neutrinos. This can be seen after some algebra. Eq.(P^) 
is now replaced by 



MiM} 



( A' + e, 
B 

\C + e2 



B 

-C + t2 



C + e2 \ 
-C + e2 
D I 



(30) 
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with A' = A + m , ei = 2mim and 62 = mm4. The above matrix can be 
diagonahzed by the following if one assumes 



tan 6* 



ei sm I 



€2 cos ( 
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< 1 



where 



tan 2(f)' 
tan 26*12 



= i?12(^12)i?23(0')^12(7r/4) 

2V2C' 



(31) 



D- A' + B'' 



2(ei cos (p' — \f2e2 sin 0') 



(32) 



5(1 + cos2 0') + A sin^ 0' + ^f2C' sin 20' - D sin^ : 
The KM matrix following from eqs.(|2^j3l|) is given by 

K{Mx) = Rl2{ei2)R23{023)Rl2{7r/A) (33) 

^23 = 0' + ^ in the above equation. The mixing angle 612 introduced by 
the symmetry breaking parameter rh has two important effects. It causes a 
departure from the exact bi-maximal mixing obtained in the symmetric limit 
and it leads to splitting among the degenerate neutrinos radiatively. Both 
these features go in the right direction in solving the solar neutrino problem. 

In order to evaluate this splitting using eq.(|), one needs to express the 
neutrino mass matrix at Mx in the charged lepton mass basis. This can be 
done using eqs.(Bl), 



/ — C23 sin 29 



X 



m 



12 



V 



-C23 COS 26112 
S23C12 



-C23 cos 26*12 
C23 sin 2^12 

-S23-S12 



■S23C12 \ 



-S23S12 





(34) 



where C23 = cos(0' + 9) ; S23 = sin(0' + 9). In the absence of the symmetry 
breaking term m, ^12 = and the above matrix has a Lg — L^ — L,- symmetry. 
Its evolution using eq.(|^) does not lead to any sphtting. With non-zero 6*12, 
one finds using eq.(H) 



-45^m^C23 sin 26*12 



(35) 
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Where 5^ is radiative correction defined in eq.(0). We have neglected the 
electron Yukawa couplings in the above equation. Note that due to the spe- 
cific texture of eq. , the splitting is determined by the /i Yukawa couplings 
rather than the more dominant tau coupling. With = m?, one obtains, 

As 2.7 10-7 eV^cas sin 2^12 

a; 

The evolution from Mx to the low scale also affects the leptonic mixing. 
This can be determined after evolving eq.(p4D to the low scale. One finds 
that if the leading terms corresponding to only r Yukawa couplings are kept 
then the mixing matrix at low scale /i remains formally the same as eq. (|33| ) : 

= i?12(^12)i?23(^23(/^))i?12(7r/4) (37) 

where 

tan6l23(/u) = tan 5123(1 + ^r) • 

It follows that the RG evolution does not change the mixing matrix appre- 
ciably compared to its form eq.(|5BD at a high scale Mx- 



Phenomenological implications of eqs. (|36| , |37D depend upon the value of 
9i2 which is determined by the symmetry breaking parameter m. We take 
it as a free parameter and consider two interesting extremes corresponding 
to very small and large 612 respectively. For small ^12, the leptonic mixing 
matrix eq. ( pT] ) has a nearly bi-maximal form. The solar mass scale in eq. ( p^ ) 
could lie in the range corresponding to vacuum oscillation solution. This 
happens for moderately large value of tan/3 which is determined by 9i2. 
More specifically, one needs ^^^^2^^ ~ 10^^ in order to obtain the solar scale 
around A5 ^10"" eV^ 

Recent observations |ll| of neutrino energy spectrum and the day night 
asymmetry are found to be less favourable for the vacuum and small angle 
MSW (SAMSW) solutions respectively. In contrast, the LAMSW solutions 
including the one with the low A5 are found to be allowed. One can obtain 
this LOW MSW solution for somewhat larger value of 612 in our case. The 
mixing matrix (^) departs from purely bi maximal mixing in this case: 



I Ci2-;il2C^3 Ci2+S12C^3 „ „M \ 

' — TT^ — TT^ ^12^23 \ 

-Sl2-Ci2C^3 -Si2+Ci2C^3 u 
TT^ ^12523 

yf2 V2 ^23 



(38) 
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where c^g = cos 6*23 (/u) ; s^g = sin 6*23 

The above matrix has the correct form to simultaneously solve the solar 
and atmospheric neutrino anomalies without violating the constraint from 
CHOOZ. The latter constraint does not allow Ou to be very large as it needs 

KM = S1243 < 0.2 

The data on atmospheric neutrinos require 

sin^ = 47^^3(1 - Kl,) = As^Jcl,{l - ^^a') ~ 0.84 - 1 

The MSW LOW solution is obtained for 

As ^ (7-20) X 10-^eV^ 
sin2 2^5 = 4X2^X^22 = (c?2 - s?2C23^)^ ~ 0.68 - 0.98. (39) 

Ranges in values of ^i2j^23(Ai) and tan/? exist which satisfy the last three 
equations simultaneously and lead to the solar scale in the range required for 
the MSW LOW solutions, e.g. 

S12 ~ 0.28; S23 — C23 ~ —/=, tan/3 ~ 20 

v2 

satisfy the CHOOZ constraint and imply 

sin^ 208 ~ 0.77; sin^ 29 a ~ 0.99; Ag ~ 1.6 x 10"^ eV^ 

when Aa ~ 4 x 10~^eV^. Since the mass splitting in the model is governed 
by the /i Yukawa coupling, it is relatively small and one needs to have large 
tan P in order to obtain the scale relevant for the LOW solution. 

4 Summary 

Solution to the neutrino anomalies, specifically the solar neutrino deficits, 
may require scales as small as A5 ~ 10~^ — 10~-'^^eV^. It is interesting to 
suppose that such a scale is generated radiatively. The most economical 
possibility in this context would be to assume that the standard model inter- 
actions themselves are responsible for generating such a small mass difference. 
This requires specific textures for the neutrino mass matrix at high scale. We 
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have discussed in this paper general conditions which ensure degeneracy of 
the masses even after radiative evolution at a low scale. It requires specific 
form, eq.(]T7|) and massless or decoupled third neutrino. 

We also argued that pseudo Dirac structure originally discussed by Wolfen- 
stein can arise from a broken discrete symmetry and presented a specific ex- 
ample based on discrete symmetry in the context of three generations. This 
example provides a nice realization of the phenomenologically successful large 
angle solutions to the neutrino anomalies. 

Acknowledgements: We thank A. Dighe and S. Lola for useful comments. 
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